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Dedicated to Pierre van Moerbeke on his sixtieth birthday. 



Abstract. We prove a number of results concerning the large A'^ asymptotics of the free energy of 
a random matrix model with a polynomial potential V{z). Our approach is based on a deformation 
TtV{z) oi V{z) to z , < t < oo and on the use of the underlying integrable structures of the 
matrix model. The main results include (1) the existence of a full asymptotic expansion in powers 
of A*'"^ of the recurrence coefficients of the related orthogonal polynomials, for a one-cut regular V; 
(2) the existence of a full asymptotic expansion in powers of A''" of the free energy, for a V, which 
admits a one-cut regular deformation TtV; (3) the analyticity of the coefficients of the asymptotic 
expansions of the recurrence coefficients and the free energy, with respect to the coefficients of V; 
(4) the one-sided analyticity of the recurrent coefficients and the free energy for a one-cut singular 
V; (5) the double scaling asymptotics of the free energy for a singular quartic polynomial V . 



1. Introduction 

The central object of our analysis is the partition function of a random matrix model, 

/OO /* OO 



^<j<k<N 

N\ n /i„, (1.1) 



7V-1 



n=0 

where V{z) is a polynomial, 

2d 

V{z) = Y,VjZ^, V2d>Q. (1.2) 
i=i 

and hn are the normalization constants of the orthogonal polynomials on the line with respect to 
the weight e'^^^^), 

P„(z)P„(z)e-^^(^)dz = Mnm; Pn{z) = Z^ + ... (1.3) 

) 

In this work we are interested in the asymptotic expansion of the free energy, 

Fn = -^\^Zn, (1.4) 

as A'^ ^ cxD. Our approach is based on the deformation rj of V{z) to z^, 

n: V{z) ^ {l-r^)z'^ + V{t-^''^z), l<t<oo, (1.5) 
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SO that 

nV{z) = V{z), T^V{z) = (1.6) 
and our main reasults are the following: 

(1) under the assumption that V is one-cut regular (for definitions see Section ^ below), we 
obtain a full asymptotic expansion of the recurrence coefficients 7„, of orthogonal polyno- 
mials in powers of A^^^, and we show the analyticity of the coefficients of these asymptotic 
expansions with respect to the coefficients Vk, k = \, . . . , 2d] 

(2) under the assumption that TtV is one-cut regular for t > 1, we prove the full asymptotic 
expansion of F^r in powers of A^^^, and we show the analyticity of the coefficients of the 
asymptotic expansion with respect to Ufc, A; = 1, . . . , 2(i; 

(3) under the assumptions that (i) V is singular, (ii) the equilibrium measure of V is non- 
degenerate at the end-points, and (iii) TtV is one-cut regular for t > 1, we prove that 
the coefficients of the asymptotic expansion of the free energy of TtV for t > 1 can be 
analytically continued to t = 1; 

(4) for the singular quartic polynomial, V{z) = (-z^/4) — z^ , we obtain the double scaling 
asymptotics of the free energy of TtV{z) where t — 1 is of the order of N~'^^'^; we prove that 
this asymptotics is a sum of a regular term, coming as a limit of the asymptotic expansion 
for t > 1, and a singular term, which has the form of the logarithm of the Tracy- Widom 
distribution function. 

In result 0, the existence of a full asymptotic expansion of the free energy in powers of N^'^ was 
first proved by Ercolani and McLaughlin |EM| . under the assumption that the coefficients of V 
are small. It was used in |EMj to make rigorous the Bessis-Itzykson-Zuber topological expansion 
|BIZj related to counting Feynman graphs on Riemannian surfaces. The approach of Ercolani and 
McLaughlin is based on an asymptotic analysis of the solution of the Riemann-Hilbert problem, and 
it is very different from our approach, which is based on the deformation equations. Also in result 
((21), our proof of the analyticity of the coefficients of the asymptotic expansion with respect to Vk 
uses an important result of Kuijlaars and McLaughlin |KMj . that the Jacobian of the map of the 
end-points of the equlibrium measure to the basic set of integrals (see Section |1] below) is nonzero. 
In result ©, the existence of an analytic continuation of the free energy to the critical point (the 
one-sided analyticity) from the one-cut side was proved by Bleher and Eynard for a nonsymmetric 
singular quartic polynomial, see the paper |BEj . where, in fact, the one-sided analyticity was proved 
from the both sides, one-cut and two-cut, and a phase transition of the third order was shown. Thus, 
result © gives an extension of the result of [BE' to a general singular V from the one-cut side. 
Observe that the analytic behavior of the free energy from the multi-cut side can be different for 
different singular V and it requires a special investigation. In result (HJ, to derive and to prove 
the double scaling asymptotics of the free energy we use and slightly extend the double scaling 
asymptotics of the recurrent coefficients, obtained in our paper |BI2j . In addition, we develop the 
Riemann-Hilbert approach of |DKMVZ] for the case when t = 1 + cN~^'^^^^~^^ , where c, e > 0. In 
this case the lenses thickness vanishes as A^~^/^ but this is enough to estimate the jump on the 
lenses by e~*^^' and to apply the methods of |DKMVZj . 

The set up of the rest of the paper is the following. In Section [21 we derive formulas which 
describe the deformation of the recurrence coefficients and the free energy for a finite A^, under 
deformations of V. Here, we make use of the integrability of the matrix model, and we refer the 
reader to excellent recent surveys of van Moerbeke |vMolj . |vMo2j on different modern aspects as 
well as the history of the matter. In Section |21 we use the deformation TtV{z) to obtain an integral 
representation of the free energy for a finite A^. In Section^ we obtain different results concerning 
the analyticity of the equilibrium measure for the g-cut regular case. In Section |21 we obtain one of 
our main results about the asymptotic expansion of recurrence coefficients in the one-cut regular 
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case. This is applied then in Section El to obtain the asymptotic expansion of the free energy, 
assuming that TtV is one-cut regular for t G [l,oo). In Section [3 we derive an exact formula for 
the limiting free energy in the case when V is an even one-cut regular polynomial. In Section |H1 we 
obtain a number of results concerning the one-sided analyticity for singular V. Finally, in Section 
iniwe obtain the double scaling asymptotics of the free energy for the singular quartic polynomial 
V. 

2. Deformation Equations for Recurrence Coefficients and Partition Function 
Define the psi-functions as 



Then 



1 V(z) 
lpn{z) = —=Pn{z)e 2 , 



The psi-functions satisfy the three term recurrence relation, 

Zipniz) = 7n+l'i/'n+l(^) + Mn{z) + 7n'(/'n-l(z) 

where 



In 



h 



n—1 



Set 



/Mz)\ 

il)i{z 

'4>2{Z 

V ■■ ) 



(2.1) 
(2.2) 
(2.3) 
(2.4) 

(2.5) 



Then ()2.3() can be written in the matrix form as 

//3o 71 
71 /?i 72 



z^{z) = (z), Q 






72 h 73 
73 /33 74 
74 /34 



V 



(2.6) 



Observe that Z^r, 7^, /3„ are functions of the coefficients vi, . . . ,f2d of the polynomial V(z). 
We will be interested in exact expressions for the derivatives of Z^r, 7^, /?„ with respect to ffe. 
Set 

Vk = Nvk, k = l,...,2d. 



Proposition 2.1. We have the following relations: 



dvk 
d-yn 



dvk 

dvk 



( [Q'']n-l,n-l 



n,n—l Tn+1 n+l,n) 



(2.7) 

(2.8) 
(2.9) 
(2.10) 



where [Q'^ 



denotes the nm-th element of the matrix ; n, m = 0, 1, 2, . . . . 
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Proof. Formula (|2.8j) is proven in Eyn . By (|2.4)) . it implies (|'2.9)) . Let us prove (|2.1()j) . Introduce 
the vector function 

- (t(if ) ^ 

As shown in |Eyn| (see also |BEHj ). it satisfies the deformation equation 



dvk 



Uk{z;n)^n{z), (2.12) 



where 



f[Q{z;k - l)]n,n-l -[Q{z;k - l)]n-l,n-l 

+ [Qiz;k-l)U -[Qiz;k-l)U,n-i 



(2.13) 



and 



From (1231), 



k-l 

Q{z;k-l) = ^z^Q''-'-^. (2.14) 

3=0 

= ^C/(^; n)^„(z), C/(z; n) = f ° ^"+3 ) . (2.15) 
The compatibility condition of (|2.13)) and (|2.15|) is 

^^i!'"^ = C/fc(z;n+ l)C/(z;n) - C/(z; n)C/fc(z; n) + J—^:^U{z;n). (2.16) 

Ot;*: 7n+l OVk 

By restricting this equation to the element 22, we obtain (|2.1()|) . Proposition 12. II is proved. □ 

We will be especially interested in the derivatives with respect to V2- For k = 2, Proposition 12.11 
gives that 

din hn 



dv2 



-il-f^l-rn+i. (2.17) 



d^n In 



^(7Li + /3Li-7'+i-/3'), (2.18) 



dv2 

= llPn-l + lll3n - ll+lPn - 7^+l/3n+l- (2-19) 



dv2 

Observe that all these expressions are local in n, so that they depend only on the recurrent coef- 
ficients with indices which differ from n by a fixed number. Our next step will be to get a local 
expression for the second derivative of Z„. 

Proposition 2.2. We have the following relation: 
52 In Zn 



ll (7I/-1 + 7^+1 + + WnHn-i + Hl-i) ■ (2.20) 



Proof. For the sake of brevity we denote (') = From ()2.17|) - (|2.19() we obtain that 

(ln/l„)" = -27^7^ - 2(5n0^ - 2-fn+lin+l = -7n (7n-l + /?n-l " 7n+l " 

- 2/?„ {ll[in-l + llfin - ll+lPn - jl+iPn+l) - 7n+l (^n + - 7n+2 " (2-21) 
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where 

In = 7n (tLi + 7^+1 +Pn + "^PuPn-l + PLi) ■ (2-22) 

From and (|2.21|) we obtain now the telescopic sum, 

N-l N-1 

(In Zn)" = J2 ^-)" = (^"+1 - = - ^0- (2-23) 

n=0 n=0 

Observe that Iq = 0, because 70 = 0, hence ()2.20|) follows. □ 
Remark: When k = 1, Proposition 12.11 gives that 

diri K r. din In, a p.. 2 2 

^^ = -/3n, 5- = Y - /?n) , = 7. " 7„+l, (2-24) 

hence 

9^ In Zn o 

-^ = 7i (2.25) 
Similar formulae can be derived also for > 3, but they become complicated. 

Remark: For the case of even potentials, equations (|2.11l - [?.l())) . as well as the statement of Propo- 
sition [23 were obtained in |FIKj . It also worth noticing that in the even case, differential-difference 
equation (|2.9)) is the well-known Volterra hierarchy whose integrability was first established in 1974 
- 75 in the pioneering works of Flaschka F , Kac and van Moerbeke |KvMj . and Manakov 'Man' , 
and whose particular case (|2.18|) is the classical Kac-van Moerbeke discrete version of the KdV 
equation |KvMj . 

Remark: For the case of the even quartic potential V{z) = V2z'^ + V4^z^, Proposition 12.21 was 
proven in [FIKOj . 

3. Free Energy for a Finite 
In terms of V2 formula (|2.2())) reduces to the following: 



d'^FN _ ^ 2 /'„2 I „2 , a2 



ll (7^-1 + 7^+1 + f3l + 2PnPn-i + Pn-i) , (3.1) 



^2 

where Fjy is the free energy, see (|1.4j) . 

The main problem we will be interested in is an asymptotics of the free energy as — > cxd. Our 
approach will be based on a deformation of the polynomial V{z) to the quadratic polynomial z'^. 
To that end we set 

W{z) = V{z) - z^, (3.2) 
and we define a one-parameter family of polynomials, 

V{z;t) = z"^ + W (^-^y, t>l. (3.3) 

Then obviously, 

V{z-l)=V{z), V{z;oo)=z^. 

It is convenient to introduce the operator r^, see ()1.5() . Then V{z;t) = TtV{z). The operators rt 
satisfy the group property. 

Proposition 3.1. 

TtTs = Tts- (3.4) 
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Proof. We have that 

Tt{T,{V{z))) = Tt{{l - S-^)Z^ + V{s-^/^z)) = (1 - t'^z'^) + (1 - S~^)t-h^ + V{t-^/h-^/^z) 

= (1 - t-h-')z^ + V{t-^/h~^/^z) = TtsiViz)). 

(3.5) 

Proposition 13.11 is proved. □ 

Let Z]\f = Z]\f{t) be the partition function (jl.ip for the polynomial V{z;t) and F/v = Fj^it) the 
corresponding free energy. 

Proposition 3.2. 

= {^Nit) [7lr-i{t)+7'N+At) + P%{t) + 2PN{t)f3N~i{t) + f3%_,{t)] - i| , (3.6) 
where ^nit),l3n{t) O'TG the recurrence coefficients of orthogonal polynomials with respect to the weight 

^-NV{z;t) _ 

Proof. By the change of variables Zj = \/t uj , we obtain from that 

ZM{t)=t^'/^ZN{t), (3.7) 

where 

/oo poo 
... Yl (uj- Ukfe~ ^J=i ^^("^^*)dni ...duN (3.8) 

oo J-oo i<j<fc<Ar 

is the partition function for 

V{u;t) =V{Viu;t) =tu^ + W{u). (3.9) 

Hence 

Fr,{t) ^ InZ^it) = H-t^y^Z^it)] = F^(t) + . (3.10) 

By (jSU), 

= -llit) [ll^iit) + fN+iit) + PUt) + 2Mt)^N-i{t) + P%^^{t)) , (3.11) 
where Jnit), I3n{t) ai^e the recurrence coefficients of the orthogonal polynomials 

Pn{u-t)=t-'''^Pn{Vtu-t) (3.12) 

with respect to the weight e"^^^"'*-* . Since 

{Vt U) PniVt U; t) = Pn+l {Vt u; t) + (3n{t)PniVi U] t) + -fl{t)Pn-l (Vt U] t) , (3.13) 

we obtain that 

ln{t) = t-^/^Jn{t) , /3„(t) = t-^/^pn{t) , (3.14) 

hence from (|3.11|) . (|3.1U|) we obtain that 

^(i^7v(t) + ^) =-^[7^_i(t) + 7^+i(i)+/3^(i) + 2/3^(t)/3^-i(i)+/3^_i(t)] , (3.15) 

which implies (|3.6|) . Proposition 13.21 is proven. □ 

We would like to integrate formula 1)3. 6|) . To that end we need an asymptotic behavior of the 
recurrence coefficients 7n,/?n, n = N — 1, N, A^ + 1 ast^oo. For a finite N it is easy. 
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Proposition 3.3. Assume that n and N are fixed. Then, as t ^ oo, 



7n(t) = ^^ + or ^/'), Pn = Oit-'/'). (3.16) 

Proof. When t = oo, V{z;t) = z^, hence 7^(00), /?„(oo) are recurrence coefficients for the Hermite 
polynomials, 

7n(oo) = y^, /3„(oo) = 0. (3.17) 

When t is finite, the orthogonal polynomials and recurrence coefficients can be obtained through 
the Gramm-Schmidt orthogonalization algorithm. Since for any moment we have the relation, 



00 , „ „ ^ /'OO 

Z' 

-00 J —00 



k^.N{z^+W{z/Vt))^^^ zV^-'dz + 0(t-i/2), (3.18) 



(pnHi) follows. □ 

From Propositions 13.21 and 13.31 we obtain the following formula for Fjs[. 
Theorem 3.4. 

Jt 



+I31,{t) + 2(3n{t)I3n-i{t) + /?|,_i(r)] -\\dT, 



(3.19) 



where 7n (''"); Pn{T), n = N — 1, N, A^+1, are the recurrence coefficients for orthogonal polynomials 
with respect to the weight e~^^^^''^\ and 

is the free energy of the Gaussian ensemble. 
Proof. From Proposition 13.31 we obtain that 

l%ir) [7^_i(r) + 7^+1 (r) + /3^(r) + 2/3^(r)/3;v-i(r) + /3^_i(t)] - i = 0{r-'/') , (3.21) 

hence the integral in (|3.19jl converges. In additon, i^^r (00) is the free energy for the Gaussian 
ensemble. Since for the Gaussian ensemble. 



hn - /ioTi . . . 7n - ^ (2iV) • • • (2iV) " ^ {2NY ' ^^'^^^ 



we obtain from l|l.lj) that 



Af-l 

'IT n\ 



Zn{oo) = N\W =P^^\\n\, (3.23) 

„=o [\/iV (2iV)"J (2Ar)A^V2 11 ^ ' 



so that F/v(oo) = F^'^"'^'^. Denote the function on the right in 1)3. 19() by F^it). Then by Proposition 
^^1^ = W + ^N+iit) + /3^(t) + 2/3^(t)/3;v-i(t) + P'n-M - \ 



d^FN{t) 



(3.24) 
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hence Fjv(t) — -FAr(t) = at + b. Since F/v(oo) = Fn{oo) = Fjy^"'^'', we obtain that a = b = 0, hence 
FN^t) = Fiy{t). Theorem 13.41 is proven. □ 

4. AnALYTICITY of THE EQUILIBRIUM MEASURE FOR A REGULAR V 

This section is auxiliary. We prove in this section a general theorem on the analyticity of the 
equilibrium measure with respect to perturbations of a regular V. The proof will follow directly 
from a result of Kuijlaars and McLaughlin |KMj . that the Jacobian of the map of the end-points of 
the equilibrium measure to the integrals {TjjNj^} is nonzero. Let us introduce the main definitions. 
We will assume that V{z) is a real analytic function satisfying the growth condition, 

lim XM^ = oo. (4.1) 

|x|^oo log \x\ 

The weighted energy of a Borel probability measure fi on the line is 

Iy(„) = -ll \og\x -y\dv{x)dv{y)+ I V{x)du{x). (4.2) 

There exists a unique equilibrium probability measure feqj which minimizes the functional Iv{i^), 

Iviueq) = min {/y(i/) : 1/ > 0, / du = 1}. (4.3) 

As shown in |DKMj . the equilibrium measure is absolutely continuous and it is supported by a 
finite number of intervals, 

SUpp I^eq = [oi , bi] . (4.4) 
The density of Ugq on the support is given by the formula 

p{x) = -l-h{x)y^R+{x), X eul^-^[ai,bi], (4.5) 
Zm 

where the function h is real analytic and 

q 

R{z) = l[[{z-ai)iz-b,)]. (4.6) 

For ^y R{z) the principal sheet is taken, with cuts on U^^-^[aj, and ^y R+{x) means the value 
on the upper cut. The equilibrium measure, = feq, satisfies the following variational conditions: 
there exists a real constant I = ly such that 

Lv{x)-^^ = l, xGU^^Jai,^,], (4.7) 



Liy{x) - < X G Ml \ Joi, h], (4., 



where 



Set 



Lu{z)= / \og\z - x\dv{x), zgC. (4.9) 

Then 

aj{z) = lv'{z)-h{z)y^), (4.11) 
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see, e.g., |l)KMV/j . This implies that 

Hz) = ^<f ^^'^ ds, (4.12) 

where F is a positively oriented contour in around {z} U?^-^ [oj, bi]. Also, for j = 0, . . . ,q, 

Tj = (h dz = 2d in. (4.13) 

' 27Ti Tr /R(i) ^ ' 



where F is a positively oriented contour in Q. around uf^;^[aj, 
A real analytic V is called regular if 

(1) Inequality is strict for all a; G \ U^^Ja^, 6^], 

(2) h{x) > for all x £ uf^iiai, bi]. 

Otherwise V is called singular. We formulate now the main result of this section. 

Theorem 4.1. Suppose V{z;t), t £ [—to, to], Iq > 0, is a one-parameter family of real analytic 
functions such that 

(a) there exists a domain 17 C C such that M C 17 and such that V{z; t) is analytic on 0, x 
[-to, to], 

(b) V{x, t) satisfies the uniform growth condition, 

mm{V{x;t) : |t| < to} / a a\ 

hm \ , ,' — = oo, (4.14) 

|x|^oo log|x| 

(c) V{z; 0) is regular. 

Then there exists ti > such that if t £ [—ti,ti], then 

(1) V{z;t) is regular, 

(2) the number q of the intervals of the support of the equilibrium measure of V{z; t) is inde- 
pendent oft, and 

(3) the end-points of the support intervals, ai{t),bi{t), i = 1,. . . ,q, are real analytic functions 
on [-ti,ti]. 

Proof. The regularity oiV{z;t) and t-independence of q are proved in |KMj . To prove the analyticity 
consider the system of equations on {oj, bi, i = 1, . . . , q}, 

Tj=26kq, j = 0,l,...,q; Nk = 0, k=l,...,q-l, (4.15) 

where Tj is defined in (|4.13j) and 

1 



Nk = — : f h{z)yjR{z) dz, (4.16) 



2m JT^ 

where Ffc is a positively oriented contour around [6^,0^+1], which lies in a small neighborhood of 
Cfc+i], so that Ffc C 17 and it does not contain the other end-points. In (|4.1(i|) it is assumed that 
the function yR{z) is defined in such a way that it has a cut on [6^,0^+1]. As shown in |KMj . 
the Jacobian of the map {[aj,6j]} — > {Tj,Nk} is nonzero. The functions Tj,Nk are analytic with 
respect to ai,bi and t. By the implicit function theorem, this implies the analyticity of ai{t),bi{t). 
Theorem 14. II is proved. □ 

When applied to a polynomial V, Theorem 14. II gives the following result. 

Corollary 4.2. Suppose V{z) = viz + - ■ ■ + V2dz'^'^ , V2d > 0, is q-cut regular. Then for any p < 2d 
there exists ti > such that for any t £ [—ti,ti], 
(1) V{z; t) = V{z) + tzP is q-cut regular. 
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(2) the end-points of the support intervals, ai{t),bi{t), i = 1,. . . ,q, are real analytic functions 
on [-ti,ti]. 

Theorem 14. II can be applied to prove the analyticity of the (A^ = oo)-free energy, 

F= lim -^In Zn. (4.17) 

If V is real analytic satisfying growth condition (|4.1|) . then the limit on the right exists, see [Joh], 
and 

F = /y(l/eq). (4.18) 

Theorem 4.3. Under the conditions of Theorem \41\ the free energy F = F{t) is analytic on 

Proof. The density of the equilibrium measure has form (|4.5() . where h{x) is a real analytic function, 
which is found by formula 1)4. 12() . By Theorem 14.11 the end-points of the support of i^eq depend 
analytically on t, hence (|4.12j) implies that h depends analytically on t, and, therefore, u^g depends 
analytically on t. Formula (|4.17|) implies the analyticity of F. Theorem 14.31 is proved. □ 

Theorem 14.31 implies that the critical points of the random matrix model, the points of nonana- 
lyticity of the free energy, are at singular V only. 

5. Asymptotic Expansion of the Recurrence Coefficients for a One-Cut Regular 

Polynomial V 

In this section we will assume that V{z) is a polynomial, which possesses a one-cut regular 
equilibrium measure. The equilibrium measure is one-cut means that its support consists of one 
interval [a,b], and if it is one-cut regular then 



di^cqix) = 7rH^)V (b - - a), x G [a,b], (5.1) 
zvr 

where h{x) is a polynomial such that h(x) > for all real x (see the work of Deift, Kriecherbauer 
and McLaughlin |DKMj ) . For the sake of brevity, we will say that V{x) is one-cut regular if its 
equilibrium measure is one-cut regular. 

As shown by Kuijlaars and McLaughlin |KMj . if V{x) is one-cut regular then there exists e > 
such that for any s in the interval 1 — e<s<l-|-e, the polynomial s~^V{x) is one-cut regular, 
and the end-points, a{s),b{s), are analytic functions of s such that a(s) is decreasing and b{s) is 
increasing. (In fact, the result of Kuijlaars and McLaughlin is much more general and it includes 
multi-cut V as well.) 

Proposition 5.1. Suppose V{x) is one-cut regular. Then there exists e > such that for all n in 
the interval 

Tl 

l-e<-<l + e, (5.2) 
the recurrence coefficients admit the uniform asymptotic representation, 

7n = 7(^)+0(iV-i), p^ = pQ)+0{N-'). (5.3) 
The functions 7(5), /3(s) are expressed as 

,„) = K£)^. ^(,) = 5(f)±M, (5.4) 
where [a{s),b{s)] is the support of the equilibrium measure for the polynomial s^^V{x). 
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Proof. For n = N the result follows from |DKMVZ] . For a general n, we can write NV = ns~^V, 
s = n/N, and the result follows from the mentioned above result from |KMj . that s~^V is one-cut 
regular, and from |DKMVZ] . The uniformity of the estimate of the error term follows from the 
result from |KMj on the analytic dependence of the equilibrium measure of s^^V on s and from 
the proof in |l)KMV/j . □ 

We can now formulate the main result of this section. 

Theorem 5.2. Suppose that V{x) is a one-cut regular polynomial. Then there exists e > such that 
for all n in the interval \5.}3\) . the recurrence coefficients admit the following uniform asymptotic 
expansion as N ^ oo in powers of N~'^ : 

k=l \ 

where f2k{s), g2k{s), k>l, are analytic functions on [1 - 

Proof. Let us remind that the proof in jPKMVZj of the asymptotic formula for the recurrence 
coefficients is based on a reduction of the Riemann-Hilbert (RH) problem for orthogonal polynomials 
to a RH problem in which all the jumps are of the order of A^~^. By iterating the reduced RH 
problem, one obtains an asymptotic expansion of the recurrence coefficients, 

oo oo 

7^ ~ 7 + 5Z ' ~ + E ^'"'dk . (5.6) 

fc=l k=l 

For a general n, let us write NV = ns~'^V , s = n/N. Then, as shown in |KMj . the equilibrium 
meassure of s^^V is one-cut regular and it depends analytically on s in the interval [1 — e, 1 + e]. 
As follows from the iterations of the reduced RH problem, the coefficients f^, Qk are expressed 
analytically in terms of the equilibrium measure and hence they analytically depend on n/N , so 
that 

oo oo 
fc=l k=l 

where fk{s), gk{s) are analytic functions on [1 — e, 1 -|- e]. We can rewrite the expansion of /3„ in 
the form 



1 \ °° / „ I i 

k 



k=l \ / k=l 




where gk{s) are analytic on [1 — e, 1 -|- e]. What we have to prove is that = g^ = 
This will be done by using the string equations. 

Recall the string equations for the recurrence coefficients, 

7n[^'(Q)]n,n-l = ^ , [V^'(g)]„n = 0. (5.9) 

where [V'{Q)]nm is the element {n,m) of the matrix V'{Q). We have, in particular, that 

[Q]n,n— 1 — 7rn [Q]nn — Pm 

[Q\,n-1 = Pn-lln + Pnjn, [Q]nn = 7^ + + 7n+l ! 

[Q\,n-1 = 7n-l7n + 7n + 7n7^+l + ^n-l7n + Pn-lPnln + P^n, 

[QXn = dn-lll + 2/3n7n + 2/3n7n+l + ^n+l7n+l + P't 

and so on. 
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Lemma 5.3. For any k> 1, the expression of [Q^]n,n-i in terms ofjj,Pj is invariant with respect 
to the change of variables 

(^0 = {7j l2n-j, Pj P2n-j-l, j = 0, 1, 2, . . .}, (5.11) 

provided n > j + k. Similarly, the expression of [Q'']n,n in terms of ■^j,f3j is invariant with respect 
to the change of variables 

(y\ = {ln+j 7n-j+l, Pn+j ^ Pn-j, j = 0, 1, 2, . . .}, (5.12) 

provided n > j + k. 

Proof. Observe that the matrix is symmetric. By the rule of multipHcation of matrices, 

[Q ]n,n—l = ^ ^ Qn.ji ■ ■ ■ Qjk-i,n—l — ^ ^ Qn,a{jk_i) ■ ■ ■ Qa{j{),n-1^ (5.13) 

where a{j) = 2n — j — 1. Observe that cr{n) = n — 1, cr(a{j)) = j and 

Qjj = Pj, Qa{j),a{j) = /?2n-i-i; Qj,j~l = 7j , Qa{j),a{j-1) = l2n-j 

This proves the invariance of [Q^]n,n-i with respect to uq. The invariance of [Q'']nn with respect 
to o"! is estabhshed similarly. Lemma 15.31 is proved. □ 

Since V'{Q) is a linear combination of powers of Q, we obtain the following corrolary of Lemma 

Corollary 5.4. The expression of ^niV {Q)]n,n-i (respectively, [V'{Q)]nn) in terms ofjj,Pj is 
invariant with respect to the change of variables Uo (respectively, ai). 

Let us 

(1) substitute asymptotic expansions (|5.8|) into equations (|5.9|1 and expand into powers series 
in N-^, 

(2) expand 7 (^) , fk (^) , P (^^) , 9k (^^) m the Taylor series at s = f , 

(3) equate coefficients at powers of N^^. 

This gives a system of equations on 7,/?, fk,gk- The zeroth order equations read 

l[V'{Qo)]n,n-l = S , [V'{Qo)]nn = 0, (5.14) 

where Qq is a constant infinite Jacobi (tridiagonal) matrix, such that 

[Qolnn = P, [Qo]n,n-l = [Qo]n-l,n =7, n G Z. (5.15) 

Equations ()5.14|1 are written as 

^(7,/3) = s, i?(7,/3) = 0, (5.16) 

where 

if WzK+l 



^Cl L 2 J / • 1 \ 

^(7,/3) = 7E^--.E/5^""-V'="H2r+i) 

j=2 k=o \ "T / 

B(.,«^|;i".E'5'-^'-v'(^,-;)(t) 



(5.17) 



j=l k=0 
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Observe that 7,/? given in (|5.4j) solve equations (|5.16|) . The A:-th order equations for A; > 1 have 
the form 

dA{^,0) dA{^,p) 

where p,q are expressed in terms of the previous coefficients, j, (3, fi, gi, . . . , fk~i, gk-i, and their 
derivatives. Here the partial derivatives on the left are evaluated at 7,/? given in (|5.4() . 

Lemma 5.5. The first order equations are 

dAirf^ , dA{j,p) 

+ = °- 

Proof. Observe that the terms with fi,gi are the only first order terms which appear at step (1) 
above. All the other terms appear at step (2), in the expansion of 7 ^^^^ , fk (^^^) ; P 

gk I — — ] in the Taylor series at s = jj. Consider any monomial on the left in the first equation 
in (EH), 

C7n+ii ■ • • In+jpf^n+h ■ ■ ■ Pn+lq ■ 

By lemma there is a partner to this term of the form 

C'Jn-ji ■ ■ ■ In-jpPn-h-l ■ ■ ■ Pn-lq-l- 

When we substitute expansions 1)5. 8|) . we obtain 

c{l{s + '^) + ..^...Us + '^,) + ..^{(3{s + ^^) + ..\...{(5{s+'^) + .^ 
and 



N' / N' / \ ' N ' / V ^ N 



N' J \" N' J y ' N ' / V ^ 

for the partner. When we expand these expressions in powers of A^~^, the first order terms cancel 
each other in the sum of the partners (in fact, all the odd terms cancel). This proves the first 
equation in (|5.19p . The second one is proved similarly. Lemma 15.51 is proved. □ 

Lemma l5 . 51 implies that /i(s) = gi{s) = for all s such that 

/ 9A(7,/3) dA{-y,f3) \ 

det I oJlU mU]^^^ (5-20) 

\ 97 a/3 / 

where all the partial derivatives are evaluated at 7(5), /3(s) given in (|5.4|) . 

Lemma 5.6. If for a given s G [1 — e,l + e], condition h5.2U\) holds, then all odd coefficients 
/2fc+i(s),92fc+i(s) are zero. 
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Proof. By Lemma 15.51 fi (s) = gi{s) = 0. If we consider terms of the order of N^"^ then we obtain 
the equations 

Indeed, the same argument as in Lemma 15.51 proves that all other terms of the third order cancel 
out. Since condition (|5.2Uj) holds, it implies that fsis) = 53(3) = 0. By continuing this argument 
we prove that all odd /2fc+i(s), 52^+1 (-s) vanish. Lemma I^TBl is proved. □ 

Lemma 5.7. Condition \5.2U\) holds for all s £ [1 — e,l + e]. 

Proof. By differentiating equations (|5.16|) in s we obtain that 



dA{-i,f3) d-f ^ dA{-f,(3) dl3 _^ 
97 ds dp ds ' 

dB{^,P) d-f ^ dB{^,(3) dp _^ 
O7 ds dp ds 

By differentiating equations 1)5. 16(1 in ti we obtain that 

dA{j,P) 07 dA{^,P) dp _ 
d-f dti dp dti 

dB{j,P) dj dB{j,P) dp 
57 dti dp dti 

By rewriting equations ()5.22|) . (|5.23() in the matrix form, we obtain that 



(5.22) 



(5.23) 



_dJ3 _ \ I ds dti 




a/3 / \ds atj 



1 



(5.24) 



Since 



det (J =-1^0, 

this implies 1)5.20(1 . Lemma 15.71 is proved. □ 

From Lemmas 15 . 61 and l5 . 71 we obtain that the odd coefficients f2k+i, 92k+i vanish. Theorem 15.21 
is proved. □ 

6. Asymptotic Expansion of the Free Energy for a One-Cut Regular V 

We have the following extension of Theorem 15.21 

Theorem 6.1. Suppose that V{z;t), t £ [—to, to], to > 0, is a one-parameter analytic family of 
polynomials of degree 2d, such that V{z; 0) is one-cut regular. Then there exist ti > and e > such 
that for all t £ [—ti,ti] and all n £ [{1 — £)N, (1 + e)N], the recurrence coefficients corresponding 
to V{z;t), admit the following uniform asymptotic expansion as N ^ 00: 

00 /_Li\°° /_Li\ 

7n~7(^;t)+E^-"/-(^5^)' /5"-/5(^;^j+E^""^-(^;^J ' (6-1) 

where 7(5; t), P{s; t), /2fe(s; t), g2k{s'-, t), k > 1, are analytic functions of s, t on [1 — e, 1+e] x [— ti, ti]. 
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Proof. Theorem 14. II implies that the equihbrium measure of V{z;t) is analytic in t S [— (see 
the proof of Theorem M.'A^ . This implies the analyticity of 7,/? in s and t. From equations 15.181 we 
obtain the analyticity of fk,gk, k > 1. By Theorem 15.21 all odd f2k+i, g2k+i vanish. This proves 
Theorem Ifi.ll □ 

Let us return to the polynomial V{z; t) = TtV{z). We will assume the following hypothesis. 
Hypothesis R. For all f > 1 the polynomial TtV{z) is one-cut regular. 

Theorem 6.2. // a polynomial V(z) satisfies Hypothesis R, then its free energy admits the asymp- 
totic expansion, 

Fn - F^""""'' ~ F + Ar-2^(2) + Ar-4^(4) + . . . ^ (6.2) 
where _pGauss defined in \S. 20\} . The leading term of the asymptotic expansion is: 



F 

where 



1-T 



2^\t) + A^\t)P\t)-\ 



dr, (6.3) 



= ffir) = -J-W, (6.4) 

and [a(r), 6(t)] is the support of the equilibrium measure for the polynomial V{z; r). The quantities 
7 = 7(r), (3 = (3{t) solve the equations, 



^(7,/3;t) = 1, S(7,/3;t)=0, (6.5) 



where 



^ ^ j = l k=Q \ I / 

,,„,,.,..(,_i),,|;^'g,,-..-V.(i -;)(-) 



k 



(6.6) 



j=l k=0 

Proof. By applying Theorem 16. ll to TtV{z), we obtain the uniform asymptotic expansions, 

7n(t)~7(^;t)+E^"'V.(^;t) , Pnit)-pl^;tUY.N-^^g,^^-tj . (6.7) 

From ()6.6|) with r = t, as t — > oo, 

^(7, P; t) = 27' + 0{t-'/^), 5(7, 13; t) = 2(3 + 0(r ^Z^), (6.8) 
hence the solutions to system (|6.5|) are 

lit) = ^ + 0(t-V2)^ ^ 0(t-V2). (g 9) 

By differentiating equations 1)6. 5() in r = t we obtain the equations, 



(6.10) 
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where p,q are expressed in terms of j{t),(3{t) and p,q = 0{t~^^'^). From this system of equations 
we obtain that ^'{t),f5'{t) = 0(1^^^"^). By differentiating equations (|H.5|) many times we obtain the 
estimates for j > 1, 

From equations H5.18() we obtain the estimates on fk,gkj 

and from the reduced RH problem, that for any K > 0, 

K 



J >0, 



(6.11) 
(6.12) 



7«(')-7(^;t)-E^'-"A(^;« 



k=l 



< C(i^)iV-2^-2rl/2 ^ 



Pnit) - 13 



n 



K 



N 



k=l 



2k 



9k 



n + ; 

N 



(6.13) 



(cf. the derivation of the estimates HA.77|) and (|A.78|) in Appendix^) Let us substitute expansions 
1)6. 7p into p. 19(1 and expand the terms on the right in the Taylor series at n/N = 1. In this way 
we obtain the asymptotic expansion, 



e^(r) ^ jlir) [7^_i(t) + 7^+i(t) + /?^(r) + 2/?^(r)/3;v-i(T) + 

oo 

~ e(T) + ^iv-^eW(r) 



(6.14) 
(6.15) 



k=l 



where 



and 



e(r) = 27^(r)+472(T)/?2(r) 



(6.16) 
(6.17) 
(6.18) 



e^^^T) <c{k)T-^/^, k>i. 

Observe that the expression (|6.14l) is invariant with respect to the transformation 

^0 = {ij 72N-j, I3j fi2N-j-l]- 

Therefore, as in the proof of Lemmas 15.51 15.61 we obtain that all odd 0(2^+^) = 0. Theorem 16. 21 is 
proved. □ 

The following parametric extension of Theorem 16. 21 is useful for applications. 

Theorem 6.3. Suppose {V{z]u), u S [— iio,wo]} is a one-parameter analytic family of one-cut 
regular polynomials of degree 2d such that the polynomial V{z; 0) satisfies Hypothesis R. Then there 
exists ui > such that the coefficients F{u), F^'^\u), F^^\u), . . . of the asymptotic expansion of the 
free energy for V{z;u) are analytic on [—ui,ui]. 

Proof. The functions are expressed in terms of integrals of finite combinations of the functions 
'y^^^l; u) , fJ^^") {I; u) , f {1; u) , g^i^ {I; u) . By Theorem 16. II these functions are analytic in u. By the 
same argument as in the proof of Theorem 16. 21 we obtain that they behave like 0(t~^/^) as t ^ oo. 
Therefore, the integrals expressing F^'^^){u) converge and define an analytic function in u. Theorem 
16.31 is proved. □ 

The following proposition is auxiliary: it gives first several terms of the asymptotic expansion of 
the free energy for the Gaussian ensemble. 
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Proposition 6.4. The constant ^pj^auss following 



expansion: 



C'(-l) + i^) 



ln2 3 IniV , , . ^^ 1 51niV 
^-^-^--J- — + (l-ln(2»))^--f2^ 

12iV3 240A^4 + 360Ar5 + [^j^e 
Proof. This is obtained from with the help of MAPLE. 



7. Exact Formula for the Free Energy for an Even V 
For an even V, /3(r) = 0, and formula 1)6.3^ simplifies, 



1 -r 



dr, 



where 



fi(r)=7'(r). 



From (|6.5j) . (|6.6|) we obtain that i? = i?(T) solves the equation 



1. 



Set 



R{t) 



Then equation (|7.3|) is rewritten as 

1 



From (f7?T|) , 
F = 



:i-t) 



1 



i=i 



dr = lim 

T^oo 



j\{l-r)arfdT + hnT 



The change of variable r = r(^) reduces the latter formula to 

f<iT) 



F = lim 



/ 2(l-r(e))CV(e)dC+-lnr-- 



or, by ()7.5|) . to 



F = lim — 

r-^oo 2 



By (El, =^(1)- By (HSI), 
(since n = N), hence 



R(T) = ^ + o(r-i/2). 



R{T) 



2T 



+ 0{T 



-3/2^ 
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When we distribute on the right in (|7.8)) . we have the following terms: 
Km 1 



'1) lim 

T^oo 



d^ + lnT 



(2) 
(3) 
(4) 



ln2-lni?(l), 

d 



<i 



2j 



^=1 



2/c 



2d 

E 

n=2 



l<:j,k<m 
j+k=n 



2j\ /2k 



-R{1) 
n 



hence ()7.8|) can be transformed into the following expression: 

d 2d 



^ In 2 Ini? X - : X - 



where R = -R(l), 



and 



j - 2 /2j 



n=2 



2 Vj 



l<j,fc<m 



Example: The quartic polynomial, 



V{z) = — + tz'^, t>-l. 



(7.10) 



(7.11) 

(7.12) 
(7.13) 



(7.14) 



The condition t > — 1 is necessary and sufficient for one cut. By H7.3(l . R = R{1) > solves the 
equation, 



hence 

By jni, (Eini), 

Thus, by ^TU} . 



R 



+ 2ti? = 1, 

-t + + 3 



Ui = t, ^2 = 0, tt;2 = 0, t«3 = t, W4^ = - . 



In2 Ini? „ „3 9„4 

F = +tR + tR^ + -R^. 

2 2 8 



(7.15) 
(7.16) 

(7.17) 

(7.18) 
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8. One-Sided Analyticity for a Singular V 

In this section we will prove general results on the one-sided analyticity for a singular V. We 
will assume the following hypothesis. 

Hypothesis S. V{z;t), t G [0,to]! is a one-parameter family of real analytic functions such that 

(a) there exists a domain C C such that M C O and such that V{z; t) is analytic on x [0, to]) 

(b) V{x,t) satisfies the uniform growth condition, 

mm{V{x;t) : < t < tp} 
lim — = oo, (8.1) 

|x|^oo log|x| 

(c) V[z; t) is one-cut regular for < t < to, 

(d) V[z; 0) is one-cut singular and h{a) ^ 0, h{h) ^ 0, where [a, b] is the support of the 
equilibrium measure for V{z]Q). 

Theorem 8.1. Suppose V{z;t) satisfies Hypothesis S. Then the end-points a{t),b{t) of the equilib- 
rium measure for V{z]t) are analytic on [0, to]. 

Proof. Set 

TJa,b;t) = —i , ^'^^'^^^ =dz, (8.2) 
' 2m Jr ^{z-a){z-b) ^ ' 

where F is a positively oriented closed contour around [a, 6] inside Vt. For t G [0, to], we have the 
following equations on a = a(t), b = 6(t): 

ro(a,6;t)=0, Ti{a,b]t) =2. (8.3) 

By differentiating (|8.2|) we obtain that at a = a(t), h = b{t), 

da 47rz Jr (z-a)/R(i) ^.m y ,/r(z) ^^Jz-a 2 

Similarly, 

dTo{a,b;t) _ h{h) dTi{a,b;t) _ ah{a) dTi{a,b;t) _ bh{b) 
db ~ da ~ 2 ' db ~ 2 ' 

Thus, the Jacobian, 

/97o OTb\ h{a)h{b){b - a) , , , 

The function To(a, 6;t) is analytic in a,b,t, hence by the implicit function theorem, a{t),b{t) are 
analytic on [0, to]. Theorem 18. II is proved. □ 

Corollary 8.2. Suppose V{z;t) satisfies Hypothesis S. Then 

(1) the function h{x]t) is analytic on x [0,to], 

(2) the free energy F{t) is analytic on [0, to], 

(3) the functions j{t),(3(t) are analytic on [0, to]. 

Proof. The analyticity of h follows from formula (|4.12l) and the one of F, from (|4.18jl . Finally, the 
analyticity of 7(t),/3(t) follows from (|5.4j) . □ 

Theorem 18.11 and Corollarv 18.21 can be extended to multi-cut V. We will say that V is g-cut if 
the support of its equilibrium measure consists of q intervals, [aj,6j], i = 1, . . . ,q. We will assume 
the following hypothesis. 

Hypothesis Sq. V{z;t), t G [0, to], is a one-parameter family of real analytic functions such that 
(a) there exists a domain C C such that R C and such that V{z; t) is analytic on O x [0, to]. 
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(b) V{x, t) satisfies the uniform growth condition, 



min{V(x;t) : < t < to} 
lim : : — ; = OO, 



|a;|— >oo 



log |x| 



(8.7) 



(c) V{z; t) is q-cnt regular for < t < 

(d) V[z; 0) is g-cut singular (with the same q as in (c)) and h{a,i) ^ 0, h{hi) 7^ 0, i = 1, . . . , 

Theorem 8.3. Suppose V{z;t) satisfies Hypothesis Sq. Then the end-points ai{t),bi{t) of the 
equilibrium measure for V{z;t) are analytic on [0,to]- 



Proof. Consider system of equations (|4.15() for V = V(z; t). As shown in |KMj . the Jacobian of the 
map 

/ : {ai,bi, i = l,...,q} ^ {Tj,Nk, j = 0,...,q] k=l,...,q-l} (8.8) 
at {ai{t),bi(t)} is equal to 



dui dbi 



X det 



1 

ax 



1 

bi 



\J Rj^{x\)dx\ . . . 

1 \ 



\jR+{x~C)dXq-i 



(8.9) 



(xi-ai) 1 (xi-61) 1 ... {xi-bq) 1 

\(Xg_l-ai)"^ (Xq„l-6l)"^ ... {Xq-l-bq)'^/ 

The determinant on the right is a mixture of a Vandermonde determinant and a Cauchy determi- 
nant. As shown in |KMj . it is equal to 

11^=1 Y\l=i{bk - flj) Wi<j<k<q(<^k - aj){bk - bj) ni<i<fc<g-i(^fc - ^j) 



i-l)''-'nU^Uixj-ak)ix,-bk) 

which is nonzero for 

ai < 61 < < a2 < • • • < bq-l < Xq-l < Oq < bq, 

and therefore has a fixed sign. Hence the multiple integral in (|8.9|) is nonzero. Now, 



dU 1 



Viz) 



dz 



1 



doi ATTi Jy {z - ai)^/R{z) 47ri \ ^/R{z) 



2uj{z) 



+ h{z) 



1 h{ai 
dz = 



z- ai 



.10) 



^.11) 



^.12) 



and a similar formula holds for Thus, the Jacobian (|8.9|) is nonzero. The functions {Tj, Nk} are 
analytic in {aj,6j},t, hence, by the implicit function theorem, {aj(t), 6j(t)} are analytic on [0,to]- 
Theorem 18.31 is proved. □ 

As a corollary of Theorem 18. 3| we obtain the following results. 

Corollary 8.4. Suppose V{z]t) satisfies Hypothesis Sq. Then 

(1) the function h{x;t) is analytic on x [0,to]; 

(2) the free energy F{t) is analytic on [0,to]; 

Proof. The analyticity of h follows from formula (|4.12|) and the one of F, from (|4.18|) . □ 
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The following extension of Theorem lH.HI will be useful for us. Suppose V{z; t) satisfies Hypothesis 
Sq. Then, as shown in |KMj . for every t G (0, to]) there exists e = e{t) > such that for any 
s £ [1 — e,e] the function s~^V{z] t) is g-cut regular and the end-points, ai{s; t), bi{s; t) are analytic 
in s G [1 — e, e]. 



Proposition 8.5. Suppose V{z;t) satisfies Hypothesis Sg. Then for any j > the functions 



d^ai{s;t) 



d^k{s-t) 



s=l 



dsi 



,q, 



i.l3) 



are analytic on [0,to]- 



Proof. By differentiating system (|4.15jl in s and setting s = 1, we obtain a linear 2q x 2q system of 
equations on 



dai{s;t) 



ds 



dbi(s:t) 



s=l 



ds 



I, 



i.U) 



The determinant of the system is calculated in H8.9() . 1)8.10(1 and it is nonzero. The coefficients of 
the system are analytic function in t S [0, to]) hence functions (|8.14|) are analytic in t G [0, to]- By 
differentiating system ()4.15|1 twice in s and setting s = 1, we obtain a linear 2q x 2q system of 
equations on 

5V(s;t) d%{s;t) 



ds"^ 



s=l 



ds^ 



l,...,q. 



i.l5) 



s=l 



The coefficients of the system are the same as for the first derivatives, but the right hand side 
changes, and it is expressed in terms of ai,bi and its first derivatives in s at s = 1, which are 
analytic in t G [0, to]- This proves the analyticity of the second derivatives, and so on. Proposition 
18.51 is proved. □ 

Let us consider next the coefficients, 'y{s;t), f3{s;t), f2k{s;t), g2k{s',t), k > 1, of asymptotic ex- 
pansions (|5.5(1 of the recurrence coefficients for the polynomial s^^V{z;t). 

Proposition 8.6. Suppose V{z; t) is a one-parameter family of polynomials of degree 2d, which 
satisfies Hypothesis S. Then for any j > the functions 



dsj 



d^p{s;t) 



ds^ 



d'f2kis;t) 



.8 = 1 



dsj 



dsi 



A; = l,2,..., 



^.16) 



s=l 



s=l 

are analytic on [0,to]. 

To prove Proposition 18.61 we will need the following lemma. 

Lemma 8.7. Suppose V{z]t) is a one-parameter family of polynomials of degree 2d, which satis- 
fies Hypothesis S. Consider the functions A(j,l3;t), -6(7, /3;t) corresponding to V{z;t). Then the 
Jacobian, 

/ dA{j,l3;t) aA(7,/3;f) \ 
Hpt d-y dl3 \ 

V ^7 dl3 / 

evaluated at ^ = 7(t),/5 = /3(t), is analytic and nonzero on [0,to]- 



^.17) 



Proof. The analyticity follows from Theorem 18.31 Let us prove that the Jacobian is nonzero. 
Consider the two-parameter family of polynomials. 



V{z;t,ti) = V{z;t) + tiz. 



22 



PAVEL BLEHER AND ALEXANDER ITS 



Then for every t £ (0,to] there exists e = e{t) > such that V{z;t,ti) is q-cut regular for any 
ti E [— e,e]. As in Proposition 18.51 we obtain that the functions 



dai{ti;t) 



dti 



dbi{tr,t) 
ii=o ' 9h 



l,...,q, (8.18) 



ti=0 

are analytic on [0,to]- By using identity (|5.24j) . we obtain that the Jacobian (|8.17j) is nonzero. 
Lemma 18.71 is proved. □ 

Proof of Provosition [KR The analyticity of 7 and (3 follows from ()5.4() . To prove the analyticity of 



d^{s;t) 



ds 



dp{s;t) 



i.l9) 



=1 ds 

let us differentiate string equations ()5.17|1 in s and set s = 1. This gives a linear analytic int £ [0, Iq] 
system of equations, whose determinant is nonzero by Lemma [8. 71 hence functions (|8.19l) are indeed 
analytic on [0,to]- By differentiating string equations 1)5. 17() in s twice we obtain the analyticity of 
the second derivatives, and so on. 

Let prove the analyticity of /2,52- By following the proof of Lemmas 15. 5( 15. HI we obtain that 
the functions /2 , 92 also satisfy a system of linear equations with the same coefficients of partial 
derivatives of A and B and an analytic right hand side. Hence /2, 52 are analytic. By differentiating 
with respect to s the system of linear equations on /2,52 and setting s = 1 we obtain a similar 
linear system for the derivatives of /2, 52; and so on. The same argument applies to f^, and their 
derivatives, etc. Proposition 18.61 is proved. □ 

Now we can prove the one-side analyticity of the coefficients of the asymptotic expansion of the 
free energy. We will assume the following hypothesis. 

Hypothesis T. V{z) is a polynomial of degree 2d such that 

(a) TtV{z) is one-cut regular for t > 1. 

(b) V{z) is one-cut singular and h{a) ^ 0, h{b) ^ 0, where [a, b] is the support of the equilibrium 
measure for V{z). 

By Proposition VA.U if V satisfies Hypothesis T, then TtV, t > 1, satisfies Hypothesis R, hence 
by Theorem 16.21 the free energy -F/v(t) of TtV admits the asymptotic expansion, 

FNit) - F^^""'' ~ F{t) + iV-2F(2)(t) + Ar-4^(4)(i) + . . . (8.2O) 

Theorem 8.8. Suppose V{z) satisfies Hypothesis T. Then the functions F{t) and F^^''\t), k>l, 
are analytic on [l,oo). 

Proof. The analyticity of F{t) is proved in Corollary 18.21 Let us prove the analyticity of F^'^^\t), 
j > 1. To that end substitute expansions 1)5. 5|) into (|3.19() . and expand the appearing functions 
7,/3, f2ki92k in the Taylor series at n/N = 1. As a result, we obtain asymptotic expansion (|8.20() . 
so that the coefficients F^'^^\t) are expressed in terms of functions (|8.16|) . By Proposition 18.61 
functions (|8.16|) are analytic on [0,to]> hence the ones F^'^^\t) are analytic as well. Theorem 18.81 is 
proved. □ 

The asymptotics of the partition function for a singular y is a difficult question. The leading 
term is defined by the {N = C!o)-free energy F, see (|4.17j) . but the subleading terms have a nontrivial 
scaling. The behavior of the subleading terms depends on the type of the singular V. The entire 
problem includes the investigation of the scaling behavior of the partition function for a parametric 
family V{t) passing through V. This is the problem of the double scaling limit. In the next section 
we discuss the double scaling limit for a singular V of the type I in the terminology of |DKMVZ] . 
when h{z) = inside of a cut. We consider a family V{t) of even quartic polynomials passing 
through the singular polynomial V. 
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9. Double Scaling Limit of the Free Energy 

We will consider the asymptotics of the free energy near the critical point of the family TtV{z) 
generated by the singular quartic polynomial V{z) = j — z"^^ 

nV{z) ^ V{z- t) = ^z'+(l--^z\ (9.1) 

We have that V{z] 1) = V{z) = | — z^, and for t > 1 the support of the equilibrium measure 
consists of one interval, while for t < 1 it consists of two intervals. We want to analyse the 
asymptotics of the free energy F]sf{t) as ^ oo and the parameter t is confined near its critical 
value, i.e. t = 1. Specifically, we shall assume the following scaling condition, 

|(t_l)Ar2/3| ^c, 
and will introduce a scaling variable x according to the equation 

t = 1 + Ar-2/32-2/3^ _ ^9 2) 

Our aim will be to prove the following theorem 

Theorem 9.1. Let -F/v(t) be the partition function corresponding to the family V{z;t) of quartic 
potentials h9. Then, for every e > 0, 

FN{t) - F^^""'^ = F'j^'^it) + A^-2F^'^S(t) + 0(iV-7/3+^)^ (9,3) 
as N ^oo and\{t- 1)A^2/3| ^ q jj^^^^ 

Fj^^{t) = F{t) + N"^F^^\t) 
is the order N~'^ (regular at t = 1) piece of the one-cut expansion I18.2(J\) . and 



F7^(t) = -logFT^((t-l)22/3iv2/3^ 
The function Fxwix) is the Tracy- Widom distribution function defined by the formulae |TWj 

Ftw{x) = exp jy {x- y)u'^{y)dy^ , (9.4) 

where u{y) is the Hasting s-McLeod solution to the Painleve II equation 

u"iy)=yu{y) + 2uHy), (9.5) 
which is characterized by the conditions at infinity jHMj . 

lim 4£L = 1, hm ^ = 1. (9.6) 

j/^-oo / _y y^co Ai (y) 

V 2 

Proof. The proof of this theorem is based on the integral representation (|3.19|) of the free energy 
wich in the case of even potentials V{z) can be rewriten as follows 



oo 



F^{t)=F^-^'' + i-fQ^{r)dT, (9.7) 
where 

e^(t) := RN{t)iRN+i + Rn^i) - \ (9.8) 
and we have used a standard notation 

ll = Rn. (9.9) 
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Note also that for even potentials all the beta recurrence coefficients are zero. Assuming the double 
scaling substitution (|9.2j) . and making simultaniously the change of the variable of integration, 

r = 1 + iV-2/32-2/3y^ 

we can, in turn, rewrite 1)9 .Tf) as 

F^{.) = + 2-/3Ar-4/3 ^^^/.y,lysyy ^N{y)dy, (9.10) 
where, we use the notations, 

F7v(x) := FN{t)\^^-^^^-2/32-2/3^ (9.11) 

and 

QN{y) ■■= &N{r)\^=l+N-y32-2/3y (9.12) 

Our next move toward the proof of theorem 19.11 is to split the integration in (|9.1fl|) into the 
following two pieces. 

Going back in the second integral to the original variable r, we have the formula. 



+ iV-2/32-2/3y) 



where 



foo 1 + iV-2/32-2/3^ _ ^ 

+ / — 2 07v(r)(ir, (9.14) 

Jl+2-2/3Ar-i 

The main point now is that we can produce the uniform estimates for the recurrence coefficients 
Rn, and hence for the function 0^, on each of the two domains of integration. Indeed, the needed 
estimates are the extensions to the larger parameter domains of the double-scaling asymptotics ob- 
tained in |B12j (the first integral) and the one-cut asymptotics obtained (in particular) in |DKMVZ] 
(the second integral). Let us first discuss the double-scaling estimates. 
Set 

50 = ^ , = 2 - 7 > (9-15) 

so that the potential (|9.1() is written as 

F(z;t) = |/ + |z^ (9.16) 

Following |BI2j , define y as 
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Then, as shown in |BI2j . 

Rn{t) = -7^ + N-y-'cii-ir+My) + N~'/'c2v{y) + O(iV-i) , (9.18) 
f2{-K)Y^^ 1/ 1 

Cl = ^ , C2 = - 



; ' 2 V2(-'^)5o, 

as ^ cxD and as long as the values of t and n are such that y stays bounded, 

\y\ < C. (9.19) 

In l|9.18jl . u{y) is the Hastings-McLeod solution to the Painleve II equation defined in (|9.5|1 - (|9.6j) . 

and 

?;(y) =y + 2^x2(y). (9.20) 
Assume that t = 1 + N^'^^^2^'^^^y. Then, by simple calculations, we have 

^ = 1 - 2i/3iV-2/3y + 0(iV-4/3y2)^ ^^-1 ^ 2-1/3 + OiN-'/^y). (9.21) 

Therefore, 

y = i/ + 0(iV"2/3y2)^ j^^j y = y±2-i/3iV-i/3_^0(iV"2/3y2), (9.22) 
if n = A'' and n = it 1, respectively. Simultaneously, 

-l + 0(iV~^/V), (9.23) 



2^ 



Cl = 2^/3 + 0{N-'^/^y), and cs = 2-^/3 ^ o(Ar-2/3y)_ (9 24) 

Let (cf. (Umi), (Eni)) 

^n(2/) := -Rn(t)lt=i+Ar-2/32-2/3y, (9.25) 

and assume that 

\y\ < C. (9.26) 

Then, we conclude from 1)9. 18(1 - 1)9. 24(1 that, as ^ 00, the recurrence coefficients Rniy), n = 
N — 1, N, N + 1, have the following asymptotics: 

Rj^(y) = 1 - N-^/^2^/^{-l)^u{y) + Af-2/32-5/3^(y) ^ 0{N-^) , (9.27) 

RN±i{y) = RNiv) T N-^/-'2'/' u'{y) + O(A^-i) . (9.28) 

To be able to use the estimates ()9.27() - (|9.28() in the ffist integral in 1)9. 14() we need them on 
the expanding domain, i.e. we want to be able to replace the inequality (|9.26|) by the inequality 

\y\ < N'. 

Proposition 9.2. For every < e < 1/6 there exists a positive constant C = C{e) such that the 
error terms in {9.21^ - ^9.2^) . which we will denote r^iy), n = N, N + 1, N — 1, satisfy the uniform 
estimates, 

\rniy)\<CN-^+^', n = N,N + l,N-l, (9.29) 
for all N >1 and \y\ < (9.30) 
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Proof. A simple examination of the proofs of |BI2j shows that the error term in (|9.18)) can be 
specified as OiN'^y^^'^). This means that, under condition 

\y\ < < e < ^, 

6 

we have from (|9.18j) the estimate 

Rn{t) = + N-^/'ci{-ir+\(y) + N-^/'c2v{y) + 0{N-'+''/^). (9.31) 

Note that the restriction e < 1/6 is needed to ensure that the droped terms are of higher order 
that iV-2/3. Estimate (|OT|) together with (jHHH-Elll) yield the following modification of ^^TFf^ 
and (UnHl) . 

Rj^{y) = 1 - Ar-V322/3(_i)A'^(y) + Ar-2/32-5/3„(y) + 0(iV-i+2^) , (9.32) 
RN±i{y) = RN{y) T N-'/'2'/' u'{y) + 0(iV-i+2^) . (9.33) 
The error 0{N^^^'^'') is produced by the second term of 1)9. 31(1 . For instance, if n = A^, we have 

N'^/^ciu{y) = N-^/^ (2^/3 + 0(Ar-2/3y)^ J^^(y) + 0{N'^/^y^)'^ 

= N~^/\{y) + 0{N-^y^) = N-^/\{y) + 0(Af-i+2^). 

Similar arguments lead to (|9.33jl . Asymptotics (|9.32jl and (|9.33jl complete the proof of the propo- 
sition. □ 

Let us now turn to the analysis of the one-cut estimates of Rn{t) which are needed in the second 
integral in (|9.14)) . These estimates can be extracted from the general one-cut expansion (|6.13|) . In 
the case of the quartic potential (|9.1j) . the first two terms of (|6.13|) can be specified as 

R^it)=R(^;t) +N'^R(^) (IL.t)+0{r^N-^), (9.34) 



.N J \N 

n = N-l,N,N + l, N ^ oo, t = to>l, 

where the coefficient functions R{X;t) and R^'^\X;t) can be found with the help of the string 
equation 1)5. 9|1 which in the case under consideration takes the form of the single recurrence relation, 

-= i2--jRn + :^Rn{Rn+l+Rn + Rn~l). (9.35) 

We also note that the change t"^/^ in the error estimate is due to the evenness of potential 

(|9.1|) . Substituting (|9.34|1 into equation (|9.35j) we arrive to the following explicit formulae for R{X] t) 
and i?(2)(A;t). 

t 



R{X; t) = -[2-t + ^{2- ty + 3Aj , (9.36) 

R^^\X;t) = -i'-^l±^^^E^4^. (9.37) 
^ ' 8 ((2-t)2 + 3A)2 ^ ' 

We of course need an extension of the validity of the asymptotics 1)9. 34() to the large domain of the 
parameter t. 
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Proposition 9.3. For every < < 2/3 there exists a positive constant C = C{5) such that 
the error terms in \9.34-^ , which we will denote rn{t), n = N, N + 1, N — 1, satisfy the uniform 
estimates, 

\rn{t)\ <Ct~^N~'^, n = N,N + l,N-l, (9.38) 

for all N>1 and t>l + l-'^^^N'^. (9.39) 
The proof of the Proposition is given in Appendix 1X1 

We are now ready to proceed with the asymptotic evaluation of the integrals in the right hand 
side of 1)9. 14[) . We shall start with the first integral, 

Fisrt we notice that, in virture of Proposition 19.21 

e^(y) = ^ _ 2^/^N-^/\\y) + 2V3Ar-2/3y + 0{N-'+^'). (9.40) 
Therefore, the Ii can be represented as 

"AT 



where 



Assuming that 



eO,(y) := ^ - 2^I^N-^/^u\y) + 2^/^N-^/^y. (9.42) 

< e < ^ , (9.43) 

we make the error term in 1)9. 41() of order o^N^"^). 

The integral in the right hand side of (|9.4H) can be, in accordance with (|9.42)) . splited into the 
three integrals, 

^-'"N-^'^j^ ^ e%{y)dy ^ In + /12 + Ii3, (9.44) 



where 



= J. (1 + A32%3,)2 ^^(^)^^- 



and 



The integral In can be estimated, up to the terms of order A^^^^ ^s follows. 

In = / {x-y){l- 2N-^'^2-^'^y + 0{N-^'''y^))dy 
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2-4/3^-4/33 f _ _ 2Ar-2/32-2/3y)rfy + 0(iV-8/3+4^) 
2 Jx 

^^^^^-4/3 + ^!^-2 ^ ^ Ar-4/3+. ^ ^-4/3+2. 



210/3 8 27/3 210/3' 

For the second integral in H9.44() . i.e. for the integral I12, we have 

/12 = -N-^ / (x - y)(l + 0(iV-2/3y))^2(y)^y 
= / _ x)u\y)dy + 0(7V-8/3+3e) 



00 



= N-^ {y-x)u^{y)dy + 0{N'^'^+^') (9.46) 
and similarly, for the third integral, 

^13 = i:N-^ / (x - y)(l + 0{N~^/^y))ydy 



2 



-N-^ j {x-y)ydy + 0{N-^'^+^^) 



= -^N-^ - lAr-2+3e + + ^.^-8/3+4.) (g_47) 

12 6 4 ^ ^ ^ ^ 
Adding the estimates (|9.45jl . (|9.46jl . and l|9.47|) . we arrive to the following, up to the order A^~^, 
asymptotic formula for the first integral in our basic equation ()9.14p . 

^1 = -|^^"'^' + {f^y - ^)^'(y)dy + ^) + hiN,e) + 0(iV-V3+4.)^ (9 48) 

where 

hiN,e) := -^iV-^/^-^^^ + ^N-'/'^^ + liV^+Se _ ^^-2+2e, (9_49) 
Consider now the second integral in the right hand side of ()9.14p . 

1 + iV-2/32-2/3^ _ r 



f 

Jl 



QN{T)dT = h. 

+2-2/3Ar-<s 



With the help of ProDosition l9.31 we can specify the general expansion (|6.15p for our case as follows 

e^(r) = Q{t) + Ar-2e{2)(^) + 0{t-^N-^), (9.50) 

where 

@{r) = 2R\l-T)-\, (9.51) 

and 

e(2)(r) =4i2(l;r)/?(2)(i.^)+^(l.^)^^^(l.^). (9.52) 

Therefore, similar to the integral Ji, we can represent the integral I21 up to the terms of order N~'^, 
as the sum of the following three integrals, 

l2=hi + h2 + hz + 0{N~^''^), (9.53) 



ASYMPTOTICS OF THE PARTITION FUNCTION 



29 



where 



'21 



1 -r 



e(r)dr, 



and 



loo = Af^2/32-2/3^ 



/23 = A^" 



4e(T)dr, 



1 -r 



'l+2-2/3Ar-* 

In analysing each of the integrals we shall recall that 
and make use of the following elementary estimate, 



e(2)(r)dr. 



(9.54) 
(9.55) 

(9.56) 



/(r)dr = f{T)dT - f{l)s - -f'{l) - -r (1) + 0{s' 



(9.57) 



which is true under the natural conditions fulfiled in the case of each of the integrals l2k ■ Applying 
(|9.57|) to the integral I21 we obtain the asymptotic relation, 



21 



1 -r 



G(r)(ir 



1 /I -r 



2 V r2 



e(r; 



2-4/3 ^-4/3+2e 



r=l 



1 /I -r 



6 V 



-e(r 



T=l 



-8/3+4eN 



(9.58) 



Observe that 



and 



1 -T 



e(r; 



T = l 



■6(1), 



r=l 



49(1) -2e'(l). 



This, together with equations (|9.51|) and (|9.36jl allows us to evaluate the coefficients of expansion 
(|9.58|1 . Indeed we have 

e(i) = ^, e'(i) = 2, 



(9.59) 



and hence 



00 1 q 



1 

12 



210/3- 

iV-2+3e ^ 0(iV-8/3+4e)^ 

Similarly, for the integral I22 we have 

/OO 1 
^eiT)dT 



(9.60) 



-xG(l)2-^/3iV-4/3+. _ ^ ^ J_0(^^ 

2 V 



^-2/32-2/3^ / ^e(T)dT 



T=l 

2 ■ 



2-2^-2+2.^0(^-8/3+3.) 



cG(l)2-4/3iV-4/3+^ + - (26(1) - e'(l)) iV-2+2e ^ 0(Ar-8/3+3e) 
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and, taking into account (|9.59j) . 



/OO 1 
^e{T)dT 



2V3- 

The estimation of the integral I23 up to the order N~'^ is very simple - we only need to use the 
first term of (|9.57|1 : 

/OO 1 
— ^e(2)(r)dr + 0(A^-^/3+.) (9_g2) 

Adding the estimates (|9.6U|) . (|9.61|) and (|9.62j) we conclude that 

/CXD 1 _ /-CXD 1 /.OO 1 _ _ 

^e(r)dT + iV-2/32-2/3^y^ _Le(r)dr + iV-2y^ i_Ie(2)(r)dr 

-/i(iV, e) + 0(iV-8/3+4.^^ (g_g3) 

where Ii{N,e) is exactly the same collection of the epsilon-depending terms as the one which has 
appeared in formula ()9.48|) evaluating the integral Ii, and which is defined in (|9.49|) . 

Substituting estimates (|9.48j) and (|9.63j) into the basic equation (|9.14|) we obtain the following 
asymptotic representation of the free energy F]\i{x), 

/OO -1 fOO 1 

— 2-G(r)dr + 7V-2/32-2/3xy ;^e(r)dr 



210/3 

Put (cf. Km ) 



F'j^^t) = F{t) + N-^F^^Ht) 
-^@{r)dT + N-^ j i-^e(2)(T)dT, (9.65) 



and consider 



F]^^{l + 2-^l'^N-^l^x). 

It is easy to see that this object coinside with the sum I21 + 122 + -^23 (see (|9.53() ) up to the following 
formal replaicment: 

N' X. 

Therefore, we can apply H9.63() and see that 

/OO -1 /.OO -| 

i-2^G(r)dr + Ar-2/32-2/3^ j ^e{T)dT 

+Ar-2 r^-^ q(2) . w _ ^-4/3 ^ ^ 
^ t2 210/3'"^ ^24 

+0(iV-8/3) (9.66) 
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This allows us to rewrite the final equation (|9.(i4|) as 

-N-'^\ogFTw{x) + 0{N-^'-'+^'), (9.67) 
which concludes the proof of theorem 19.11 □ 

Remark. In terms of the partition function equation 1)9. 67() reads 

ZN{t) 



7Gauss 



FTw[{t-l)2'''N''')Z'^^{t) 



yi(l + 0{N-^l^+')y (9.68) 
where e is an arbitrary positive number. 

Appendix A. The proof of proposition 19.31 

Let us remind the basic steps of the Riemann-Hilbert approach to the asymptotic analysis of 
orthogonal polynomial following the scheme of |DKMVZ] . 

The principal observation ( FIKj; see also |BIlj and |DKMVZ] ) is that the orthogonal polynomials 
Pn{z) admit the representation, 

Pn{z)=Ynll{z), (A.l) 

where the 2x2 matrix function Yn{z) is the (unique) solution of the following Riemann-Hilbert 
(RH) problem. 

(1) Y{z) is analytic for z G C \ M, and it has continuous limits, Yn+{z) and l^_(z) from above 
and below the real line, 

Yn±{z) = lim Yn{z'). 

z'^z, itim z' >o 

(2) Yn{z) satisfies the jump condition on the real line, 

Yn+{z)=Yn-{z)G{z), (A.2) 

where 

/i -NV(z)\ 

G(.) = (; ^ ) (A.3) 

(3) as z ^ oo, the function Yn{z) has the following uniform asymptotics expansion: 



Yn{z) 

where 



^ oo in) \ 

^ k=l / 



0-3 



1 

In addition to equation (|A.H) . the recurrence coefficients can be also evaluated directly via 
Yn{z). In fact, we have that 

Rn = {m^i\2{mf\^, (A.5) 

where the matrix rnf'^ is the first coefficient of the asymptotic series (|A.4|) . Equation ()A.5|) reduces 
the question of the asymptotic investigation of the recurrence coefficients Rn{t) to the question 
of the asymptotic solution of the RH problem (1-3). In the case of a fixed t > 1, this analysis is 
performed in |DKMVZ] . In fact, in ^DKMVZj the asymptotics is evaluate for a generic fixed real 
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analytic potential V{z). The approach of |DKMVZ] consists of a succession of steps which, in the 
end, yields a reduced RH problem in which all the jumps are of the order A^~^ (cf. the proof of 
theorem 15 ■2j) . In the relevant for our analysis one-cut situation, these steps are described in detail 
in |EMj . In what follows we will repeat the construction of |EMj specifying its principal ingredients 
for the case of the potential (|9.1|) and showing how it can be modified in order to cover the extanded 
range of parameter t, i.e. assuming t > 1 + 2~'^^^N^^ . 
Stepl. {g-function deformation). Define 

\n{z - s)p{s)ds, (A.6) 

where [—2:0,-^0] and p{s) are the support and the density of the equilibrium measure (|5.1|) . respec- 
tively. More precisely, [— zo,2;o] and p{s) minimise the functional (|4.2j) where V{x) is replaced by 
\V{x), and A = n/N, n = N -l,N,N + 1. In the case of the potential 1)9. the point zq and the 
function p{s) are given by the equations (see, e.g. |BPSj ). 



zo = 2 Q (2 - t + V(i - 2)2 + 3A)y = 2i?i/2(^. (A.7) 



where 



and 



p{s) = ^{h^ + h2S^)^z[^^^^{b, + h2s')y7^^ , (A.8) 
60 = ^ (2t - 4 + V(i-2)2+3A) , (A.9) 

b2 = ^- (A.IO) 

The branch of y^z^ — Zq is defined on C \ [—zq, zq] and is fixed by the condition z"^ — Zq > if 
z > zq. The branch of ln(z — s) is defined on C \ (—00, s] and is fixed by the condition arg(z — s) = 
if z > s. 

Assume that t > to > 1 and denote, 

Vx{x) = jV{x). 

Then the function g{z) satisfies the following characteristic properties (cf. 1)4. 7|) - 1)4. 8(1 ) which 
underline the importance of g{z) for the asymptotic analysis of the RH problem (1-3). 

• The function g{z) is analytic for z G C \ (— oo,zo] with continuous boundary values g±{z) 
on (—00, zo]. 

• There is a constant I such that for z G [— ^Oi ^0]^ 

g+{z) + g_{z)-Vx{z) = l, (A.ll) 

and for z G M \ [— zq, zq], 

g+{z) + g^{z)-Vx{z) <l. (A.12) 

• Denote 

p(z) :=5+(z)-5_(z). (A.13) 
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Then, for z G [— zqj 2^0] , 



20 



V{z) = 2m / p{s)ds, (A.14) 



z 



and this function possesses an analytic continuation to a neighborhood of {—zq, zq). More- 
over, for every < d < zq/2 there is a positive number po such that 



d 



da 

for ah s G [—zq + d,zo — d] and t >to > 1 
For z > zo, 



Repis + ia) = ^{bo + b2S^)J zl - > po > (A.15) 
(7=0 A ' 



and for z < —zq^ 



p{z) = 0, (A.16) 
p{z) = 2m. (A.17) 



• as z — > 00, 



g{z)=lnz + o(^^y (A.18) 
We also notice that there is the following alternative representation of the function g{z), 

^^'^ = -\iy^ + b^''^\f^^ods + lvx{z) + ^. (A.19) 

Having introduced the function g{z) and the constant I, we define the first transformation, 
Y{z) ^{z) of the original RH problem, by the equation. 

In terms of the function $(2;) the RH problem (1-3) reads as follows. 

(!') $(z) is analytic for z e C \ M. 

(2') ^{z) satisfies the jump condition on the real line, 

^+{z) = ^.{z)G^{z), (A.21) 

where 

/ -np{z) n(g+{z)+g-{z)-Vx-l)\ 

= ( e-W ) ^^-''^ 

(3') as z — ^ 00, the function ^{z) has the following uniform asymptotics: 

^{z) = I + o(-), z^oo (A.23) 



(which can be extended to the whole asymptotic series). 
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C(u) 





c(i) 

Figure 1. The contour T. 
Observe that, in virtue of (|A.16|) and (|A.17|) . we have 



1 g"(5+(2)+9-{^)-^A-0 

1 



for z £R\[-zo,zo], 



and in virtue of (lA.llll 



e -np(2:) ^ 

Gq>{z) = [ ^ ) , for z G [-zo,zo]. 



(A.24) 



(A.25) 



Step2. {Second transformation <I> ^^^^) Next we introduce the lens-shaped region 
around {—zo,zo) as indicated in Figure 1 and define ^^^\z) as follows 



(i) for z outside the domain 



(ii) for z within the domain $7^") (the upper lens) 



1 oy 

_g~np{z) ^ j 1 



(A.26) 



(A.27) 



(ii) for z within the domain O^'^ (the lower lens), 



(A.28) 



(We note that the function p{z) admits the analytic continuation to the domain $7.) 

With the passing to ^^^\z), the RH problem (1' - 3') transforms to the RH problem posed on 
the contour F consisting of the real axes and the curves C^"^ and C^^^ which form the boundary of 
the domain 0,, 

n = C7(0 - 

(see Figure 1). We have. 
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where 



' fl gn{g+{z)+g.{z)-Vx-l) 

1 



1 

g-np(2) ^ 

1 

e"P(2) 1 



1 
-1 



for z G M \ [-zo,zo], 
for z £ C("), 
for z G C«, 
for z G [-^0,2^0] 



(A.29) 



(A.30) 



(3") as z ^ 00, the function ^^^\z) has the following uniform asymptotics: 



/ + 



00 



(A.31) 



(which can be extended to the whole asymptotic series). Indeed, in view of the equations (|A.26|) 
- ()A.28|) defining the function ^^^\z), the properties (!') - (3') of the function ^{z) and equation 
HA.24I) , we only need to explane the last line of equation HA.30|) . The latter is a direct consequance 
of equation (|A.25|) and the elementary algebraic identity, 



1 
-1 



1 





3"P I 



-np ^ 

n p"P 



-np ^ 



Step3. {The construction of a global aproximation to ^^^\z)) The point of the transformation 
of the original Y - RH problem (1 - 3) to the <I> - RH problem (1" - 3") is that in virtue of the 
inequalities ()A.12|) and ()A.15|) . the jump matrix (z), for z ^ ±zo, is exponentially close to the 
identity matrix on the part T \ [—zq,zq] of the jump contour F, so that one can expect that, as 

^ 00, n = iV - l,iV, + 1, and \z ± zq] > 6, 



^W(z) ^ $M(z), 
where ^^°°\z) is the solution of the following model RH problem. 

(1'") is analytic for z G C \ [-zq, zq]. 

(2'") ^>(°°)(z) satisfies the jump condition on {—zq,zq) 



(A.32) 



^^r^(z) 



<^^f°\z) 



1 
-1 



(A.33) 
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w=w(z) 





1 1 \ 








\ ^ J 


ZQ j 





2n/3 



'0 



-27C/3 



Figure 2. Decomposition of B^. 



(3"') as z ^ oo, the function ^^°°\z) has the following uniform asymptotics: 

r 



I + O 



z 



oo 



(A.34) 



(which can be extended to the convergent Laurent series at z = oo). The important fact is that 
this Riemann-Hilbert problem admits an explicit solution: 



a{z) = 





fa-1 


a— a ^ 


) = L 


-U 


2i 
a+a 




-2i 


2 


Z- Zq^ 
Z + Zqj 


r. 


a{oo) 



1. 



(A.35) 
(A.36) 



In order to prove and specify the error term in estimation ()A.32|) we need to construct the 
parametrix of the solution ^^^\z) near the end points izQ. Let i?^ denote a disc of radius d 
centered at zq, and let us introduce the change-of-the- variable function w^z) on by the formula, 



'w{z) 



2/3 



{-2giz)+Vxiz) + lf/^. 



In view of equation ()A.19|) . the function 'w{z) can be also written as, 



w{z) 



2/3 



2 f'^ 

Y / (&o + b2s' 



Zn ds 



2/3 



(A.37) 



(A.38) 



which, taking into account that |6o + > cq > for all t > 1, implies that, for sufficiently small 
d, the function w{z) is holomorphic and in fact conformal in the disc B^, 



w{z) = '^Wk{z - z^f, z e Bd 



k=l 



We shall assume that the branch of the root ( )^/^ is choosen in such a way that 

wi>co>0 for ah t>l and iV > 1. 
We also note that, for sufficiently small d, the following inequality takes place, 

l^t'l-^)! > Co, for all z G S^, t >1 and > 1, 
where Sd denote the boundary of B^, i.e the circle of radius d centered at zq. 
Let us decompose Bd into four regions (see Figure 2), 



(A.39) 

(A.40) 
(A.41) 
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where 



(A.42) 



B 



(1) 



B 



(2) 



2-K . 

z e Bd:0< arg w{z) < f ' 



27r 

z £ Bd : — < arg w{z) < vr 



B 



(3) 



z £ Bfi : —IT < aigw(z) < 



2tt 



B 



(4) 



2-K . , 

z £ Bd : — — < arg w{z) < 
o 



We shall assume that the parts of the curves C7("'') which are inside B^ coincide with the relevant 

(k) 

partys of the boundaries of the domains B^ ' . Let us also introduce the standard collection of the 
Airy functions, 

yo{z) := Ai (z), yi{z) := e'^^/^Ai (e-^^^/^^z) y2{z) := e"^/^Ai (e^^^/^^z) (A.43) 
We will now define the approximation (parametrix) <I>(^o)(z) within Bs by the following equation, 



lor z £ B 



(1) 



(2) 
d 1 



^f.{n'l^w{z))i^\ 0^gf«.^/^(.V3^ forzGS 
^«.(n2/3^(z)) (^1 0^gt-«/2(.)<x3^ forzei?f 

where the model functions ^'^■'^''(z) are the matrices, 



(A.44) 



for z £ B 



(4) 



yoiz) iyi{z) 
y'oiz) iy[{z) 



yoiz) iy2{zy 

y'oiz) iy'2{z))' 



and the gauge matrix multiplier E{z) is 



E{z) 



vr I . _i 



-a 
-ia 



wi^^z). 



(A.45) 
(A.46) 

(A.47) 



We note that, as it follows from (|A.36|) and (|A.39|) . the matrix- valued function E(z) is analytic 
in the disc Bd- 
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Figure 3. The contour Tq. 




We are now ready to define an exphcit global approximation, ^^^\z), to the solution ^^^\z) of 
the RH problem {!'" - 3"'). We take 



[ a3$(^o)(-^)fT3 for z £ i-Bd) 



(A.48) 



To see that these formulae indeed provide an approximation to the solution (^^^\z) we consider 
the matrix ratio, 



X{z) := f$(^)(z) 



-1 



(A.49) 



Due to equation HA.33|) and the definitions HA.44|) of the parametrix <I>(^'')(z), the function X{z) 
has no jumps across the interval {—zq + d,ZQ — d) and inside the discs B^ and {—Bd). It is still 
have jumps across the contour 



(-00, -ZQ -d]U i-Sd) U 



(u) 



UC^^^USdU [zo + d,+oo), 



(A.50) 



where Cq ' are the parts of the curves C*-"''^ which lie outside of the discs Bd and {—Bd). The 

curves Cq"''^ can be taken as straight lines. The contour Tq is shown in Figure 3. The matrix-valued 
function X{z) solves the following RH problem posed on the contour Tq. 

(1^) X(z) is analytic for z £ C \ Tq, and it has continuous limits, X^{z) and X^{z) from the left 
and the right of Tq. 

(2^) X{z) satisfies the jump condition on Tq 



X+{z) = X^{z)Gx{z), 



(A.51) 



ASYMPTOTICS OF THE PARTITION FUNCTION 



39 



where 



for z G M \ (— zo — (i, zo + c?)) 



for z e C7j"\ 



for zed 



(0 

' 



for z £ Sd, 
for z G i-Sd) 



(A.52) 



(3'') as z ^ oo, the function X{z) has the fohowing uniform asymptotics: 



X{z)=I + 0( - 



z ^ oo 



(A.53) 



The important feature of this RH problem is that the jump matrix Gxiz) is uniformly close to 
the identity matrix as — > oo. Indeed, using the known asymptotics of the Airy functions and 
inequality HA.41|) one can check directly that the functions ^^^°\z) and ^^°°\z) match on the circle 
Sd, and the uniform estimate, 

\Gx{z) - 1\ for all z £ SdU {-Sd), t>l, and A > 1, (A.54) 

takes place. Simultaneously, we observe that as z runs over M \ {—zq — d, zq + d), we have 

Q ^ ^n(g+(z)+g.iz)-Vx-l) = g" ^ Co +''2 ) V'^' " ^0 <^^) < g"^^"^', (A.55) 

where the positive constant cq can be choosen the same for all t > 1 and A > 1. Therefore, we 
conclude that 

\Gx{z) -I\< Ce^^^o^', for ah zeR\ {-zq -d,zo + d), t>l, and A > 1. (A.56) 

Finally, inequality HA.15|) indicates that on the segments Cq"^ and Cq^ , if they are choosen close 
enough to the real line, the estimate 



\Gx{z)-I\ < Ce 



~Nco 



for all z G C^"^ U C^'^ t>to>l, and A > 1. 



(A.57) 



holds. 

Unlike the estimates (|A.54|) and (|A.56|) . estimate (|A.57|) can not be extended to t > 1. However, 
a slightly weaker version of it is valied for t > 1 + 2-^/3 iV"^ with 6 < 2/3. To see this, let us 
analyse more carefully the behavior of the function Re p{z) near the real line. To this end let us 
notice that, in addition to (jA.lSf) we have 



—T Rep[s + la^ 
da'^ 



(7=0 



0, \fz G {-Zo,Zo), 



and hence 
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Rep{z) = a { ^Rep{s + ia) ) + O(cr^) 



a=0 

= 0(^0 + h2S^)^4 - + 0{a^), z = s + ia e C^""^ U C^^\ (A.58) 
By a straightforward calculation one can check that 

—N-^ <bo<l, Vt > 1 + 2-2/3Ar-'5. 
Therefore, equation (|A.58|) yields the estimates 

nRep{z)>coaN^-^(^l + o(^a^N^^y z = s + ia e C^''\ (A.59) 

and 

nRep{z) <coaN^~^(l + o(a^N^'^y z = s + ia £ C^^\ (A.60) 

with some positive constant cq. 
If we now choose C^"''^ so that 

\lmz\ = \a\ = iV-^/^ z e C^""^ U C^^\ (A.61) 

and assume 

then (|A.59|) and (|A.60|) would imply 

nRep(z) > cqN'^/^-^ z = s + ia e C^''\ (A.62) 

and 

nRep{z) < -coiV2/3-'5, ^ = g + g ^W. (A.63) 

(We follow the usual convention to use the same symbol for perhaps different positive constants 
whose exact value is not important to us.) These inequalities in turn yield the following modification 
of estimate (|A.57p . 

\Gx{z)-I\<Ce-^'>^'^'~\ for ah z G C^"^ U (A.64) 
t > 1 + 2~'^/''^N-\ and N>1, 
which together with HA.54|) and HA.56|) lead to the conclusion that 

(J 

\\Gx - I\\L°°(ro)^ IICx - -?'||L2(ro) < — , (A.65) 

for ah t > 1 + 2^2/^iV"^ and > 1 

This means that the needed extention of the basic uniform estimate of the jump matrix has been 
almost obtained. What is left is the control of the t - dependence of the estimate. This can be 
achieved as follows. 

Let us attach the subscript "Gauss" to all the relevant objects, i.e. the equilibrium measure, 
the model solutions, etc., which correspond to the gaussian potential, Vcauss 

(z) = z^. By the very 
nature of our approach, as t — > oo, all the main ingredients of the above scheme, i.e. 
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converge to the respective Gauss - quantities, i.e. to 

5Gauss(^), ,'J>SL(^) ^GaL(^)- 

Moreover, the fohowing inequahties for the jump matrix of the X - RH problem can be estabUshed 
by a straightforward calculations. 



I\ < 



for all z G M \ {—zq — d,xo + d), 



\Gx{z)G-^l 



z) -I < 



for all z e Cj") U C^'^ 



t > 1 + 2^2/3^-5^ ^^^^ N>1, 
t>l + 2-'^/^N-\ and iV > 1, 



(A.66) 



(A.67) 



(A.68) 



Put 



X{z) = X{z)X^ljz). (A.69) 

The function X{z) solves the RH problem on the same contour Tq as the function X{z) and with 
the jump matrix, 

G^iz)^Gxiz)G-/^_iz). 
The inequalities (|A.66l - lA.68|) yield then the following modification of estimate ()A.65|) . 



G 



(A.70) 



ll^x - Ih^iFo), \\Gx - -^llL2{ro) < ^> 

for all t > 1 + 2^2/3^-5^ jyyi 

The proof of Proposition 19.31 can be now completed in the usual way, by iterating the X - RH 
problem(cf. |DKMVZ] and |EM^ ) . Indeed, by iterative arguments, we can see that for any K > 0, 



K 

X(z)-I-^N-'h[^;t,. 



k=l 



< 



C{K) 



tNK+\l + \z\)' 



and also 



fk[^;t,z 



< 



G{k) 



(A.71) 



(A.72) 



t{i + \z\y 

Denote m^, rrai^causs and fhi the matrix coefficients of the terms l/z in the asymptotic series at 
z = oo of the functions <I>(°°)(z), -'^Gauss(-z) and X{z), respectively. Then, for the coefficient m!'^^ 



of series ()A.4|) we will have from ()A.2n|) . (|A.49|) . and ()A.69|) that 



(A.73) 
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In virtue of estimates ()A.71|) and ()A.72|) we have that 

K 



k=l 



< C{K)N-^-^t-\ 



Tk [^■,t)\<C{k)t-\ 



(A.74) 
(A.75) 



while 



-0-3 



-0-3 



m 



1, Gauss" 



(A.76) 

(n) 

Observe now that the matrices m^, ^TT-^Gauss' "^1 Gauss evaluated explicitly. Indeed, the 



first two can be obtained from HA.35|) . taking into account that zq, Gauss = v2A, and the third one 
follows from the fact that the normalizing constants /in, Gauss are known - see (|3.'22p . Therefore, 
performing the calculations indicated, we derive from equations HA.74|) . HA.73|) . and HA.5|) the 
following estimates for the recurrence coefficients 



K 



Rn{t) 



fc=l 



n 



< C{K)N- 



fk{^;t)\<c{k)t-\ 



(A.77) 
(A.78) 



for all t>l + 2-'^I^N- 



n 



N-l,N,N + l, and iV > 1, 



Finally, repeating the arguments we used in the proof of Theorem 15.21 we conclude that the odd 
coefficients in the series from ()A.77|) are actually absent. The Proposition 19.31 follows. 
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